Theory of electrostatically induced shape transitions in carbon nanotubes 
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A mechanically bistable single-walled carbon nanotube can act as a variable-shaped capacitor 
with a voltage-controlled transition between collapsed and inflated states. This external control 
parameter provides a means to tune the system so that collapsed and inflated states are degenerate, 
at which point the tube's susceptibility to diverse external stimuli- temperature, voltage, trapped 
atoms - diverges following a universal curve, yielding an exceptionally sensitive sensor or actuator 
that is characterized by a vanishing energy scale. For example, the boundary between collapsed and 
inflated states can shift hundreds of Angstroms in response to the presence or absence of a single 
gas atom in the core of the tube. Several potential nano-electromechanical devices can be based on 
this electrically tuned crossover between near-degenerate collapsed and inflated configurations. 

PACS numbers: 61.48.De, 85.85+j, 85.35.Kt, 64.70.Nd 



The equilibrium cross-sectional shape of a nanotube is 
controlled by a competition between elastic and surface 
energies. Sufficiently large-diameter tubes prefer the col- 
lapsed state, which captures the surface energy of the 
now-touching interior surfaces. Three distinct stability 
regimes can be defined in terms of the radius of the in- 
flated, cylindrical state: below R\ only the inflated state 
is stable; between Ri and i?2 the inflated state remains 
stable, but the collapsed state is metastable; above R2 
collapse is stable and inflation is only metastable 1 -"— . In 
the region of bistability above R%, transitions between 
these two configurations can propagate down the axis 
of the tube^. Here we show how the highly deformable 
conductive sp 2 sheet can act as a non-linear, variable- 
shape capacitor wherein electrostatic interactions within 
a charged tube shift i?i and R2 to favor the inflated state. 
A tube can be tuned by external voltage to a critical point 
at which the inflated and collapsed states are degenerate, 
producing a divergent susceptibility to diverse external 
stimuli and creating a regime of exceptionally sensitive 
nonlinear nano-electromechanical response. 

Consider an infinitely long single-walled tube of ra- 
dius R wherein an external voltage U (relative to infin- 
ity) imparts additional charge q onto the tube wall. We 
model the covalent interatomic interactions of this sys- 
tem with a Tersoff-Brenner potential*^, the non-bonded 
interactions with a Lennard- Jones term, and the electro- 
static contribution via a screened Coulomb interaction 
and charging energy. The parameters e and a of the 
Lennard- Jones interaction 4>lj{t) = e(^-) 12 — 2e(^) 6 
are a — 0.383 nm and e = 2.39 meV 1 -. We exclude in- 
teratomic distances r < 0.3 nm from the Lennard- Jones 
term to protect the integrity of the covalent interactions. 
For computational efficiency, we neglect non-bonded in- 
teractions for r > ro nm and shift the Lennard- Jones po- 
tential upward by a small term linear in r so that the en- 
ergy and force vanish at the cutoff distance^. In any real 
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FIG. 1: Charge distribution over the collapsed configuration 
of (20,20) tube, showing the accumulation of charge in the 
bulbs. Dark blue and light green lines correspond to U — 1 
and U = 2 V at k s = 0.1 nm -1 , while black lines correspond 
to k s = 0.05 nm" 1 at U = 1 V. 

system, the induced charge on the tube wall is screened 
by countercharges in the environment. In addition, coarse 
graining of the induced charge onto atoms requires an 
empirical on-site Coulomb self-interaction. Hence we in- 
troduce an on-atom self-energy that is linear in the in- 
duced charge and use a Thomas-Fermi screening at long 
distances: 

s (r ; ) = Aq i + £-^-£ . (1) 
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The empirical parameter A — 7.8 nm" 1 is fitted to results 
when k s = 0^, but the precise value of A has only mi- 
nor effects on our main results. Thomas-Fermi screening 
provides a reasonably accurate model of screening due to 
e.g. an electrolyte. In 10 _1 to 10 -4 M NaCl solutions 
the Debye screening length varies from 30 to 1 nm^; we 



2 



use k s = 0.2,0.1,0.05 nm l . The net electrostatic self- 
interaction is Es — \ EfcLi 9i^s( r i)- ^ fixed applied 
voltage ?7, the tube when collapsed holds less charge than 
when inflated, yielding a contribution Eq = —qU due to 
work against this external voltage. The total energy is 
then a sum of covalent (Etb), non-bonded (E^j), and 
electrostatic (E^ + Eq) contributions. The energy differ- 
ences AE quoted below are always relative to the energy 
of the corresponding inflated state at the same voltage, 
so that AE = corresponds to degeneracy. Equations 
of motion (under periodic boundary conditions) are in- 
tegrated with a Verlet algorithm, incorporating viscous 
damping to relax to static equilibria. 

A constant-voltage condition models a tube connected 
to an electrode with a quasi-ohmic contact; (a constant- 
charge condition could model a tunneling contact, with 
a distinct type of collapse/inflation transition dynamics 
described later). When the tube is collapsed, the charge 
accumulates in the bulbs, as shown in Fig. Q] To find the 
charge distribution, we write the potential 0s (r^) given 
by ([T]) for each atom 1 < i < N in the unit cell and solve 
the system of self-consistent linear equations: 

J=N n=oo e - ksPijn 
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with respect to unknown charges qi at a fixed electrostatic 
potential <fe:( r i) — U , where pij n is the distance between 
atoms i and j separated by n unit cells. The primed sum 
omits n = 0. 

Fig. Q] shows the excess charge per atom. As regions of 
high curvature, the bulbs accumulate charge, more so at 
higher voltages. More weakly screened charge distribu- 
tions require less total charge to sustain a given voltage. 
For strong screening, the on-site self-energy in Eqn. 
is more important compared to the inter-atomic interac- 
tion, so the system becomes slightly more sensitive to the 
empirical parameter A^. 

The inflated state holds more charge than the collapsed 
state at fixed voltage, therefore it has a higher capaci- 
tance and is favored under increasing charge. The upper 
panel of Fig. [5] depicts voltage-controlled shape transi- 
tions for (20,20), (30,30) and (40,40) tubes. All three 
systems are initialized to a collapsed state at U = 0; (col- 
lapse is metastable for (20,20) or (30,30) and stable for 
(40,40)). Charging not only favors the higher-capacitance 
inflated state, but also decreases the barrier against in- 
flation. At C/crit(i?) the kinetic barrier against inflation 
disappears and the tube inflates. The inset in the figure 
depicts this transition for the (40,40) tube. As the volt- 
age increases, the bulbs expand and the flattened interior 
shrinks. At [/ crit w 39 V the tube snaps open. The mid- 
dle and lower panels of Fig. [2] demonstrate influence of 
the screening parameter k s on the energy AE and charge 
difference Aq correspondingly. The effect is robust across 
a wide range of screening lengths. Of particular interest 
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FIG. 2: Upper panel: Net energy AE = AE T b + AE LJ + 
AE S + AEq per unit length of (20,20), (30,30) and (40,40) 
tubes as a function of the applied voltage U, where AE — 
corresponds to the inflated state at the same voltage. Dis- 
continuities correspond to abrupt transitions into the inflated 
state. The inset shows changes of the shape of a (40,40) tube 
under external voltage. Middle panel: AE for a (40,40) tube 
when k s = 0.2, 0.1, 0.05nm _1 . Lower panel: charge differ- 
ence between the inflated and collapsed configurations for a 
(40,40) tube; the near-linearity indicates that the end-state 
tube shapes are relatively insensitive to charge (i.e. the ca- 
pacitances of inflated and collapsed states are each roughly 
constant across this range of voltage). 



are the more modest voltages needed to tune the system 
to the degeneracy point AE = 0, a condition in which 
the system becomes exceptionally sensitive to external 
perturbations, as described below. 

Fig. [3] depicts the various contributions to AE (solid 
line), the energy difference between collapsed and inflated 
states, as functions of U for the (30,30) tube. The elastic 
contribution (blue dashed line) is always positive, since a 
uniform circular cross-section minimizes the total curva- 
ture energy. Conversely, the Lennard- Jones contribution 
(red dots) is always negative, since the collapsed state 
allows closer approach of opposing surfaces. Since the 
collapsed state has a lower capacitance, the Coulomb self- 
energy AEy, (pluses) is negative (at fixed voltage). The 
charging contribution AEq (triangles) associated with 
motion of excess charge on or off the tube is positive, 
due to the same difference in capacitance. 

For tubes with R> R 2 , like (40,40), the external volt- 
age can be tuned so that collapsed and inflated states are 
degenerate, i.e. AE = 0. Since this [/ deg cn < f/ C rit, the 
two states are separated by a barrier. If boundary condi- 
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FIG. 3: Decomposition of the energy difference AE between 
collapsed and inflated states for a charged (30,30) tube into 
covalent /elastic (Etb), Lennard- Jones (Elj), Coulomb (-Be), 
and charging (AqU) contributions, all measured per unit 
length for k s = O.lnm -1 . 



tions that eliminate this barrier can be imposed, then the 
susceptibility against perturbations that favor collapse or 
inflation will diverge at this critical voltage. For example, 
a tube with R > i?2 held inflated at one end and collapsed 
at the other (as shown in Fig. 2]) must contain a transition 
region. The kinetic barrier between degenerate configura- 
tions is thereby eliminated, since the transition region can 
move freely along the tube axis. Such pinning has already 
been produced experimentally, even in the first discovery 
paper on nanotube collapse^, since the rigid end-cap of a 
tube holds the end of the tube open even after the interior 
collapses. This pinned-open/pinned-closed configuration 
provides many possible modes of device operation; these 
can be classified as those that operate around the degen- 
eracy point versus those that operate at the degeneracy 
point. 

Devices that operate around the degeneracy point: 

By sweeping the applied voltage across f/degen, a dou- 
bly pinned nanotube can act as a electrical-mechanical 
transducer, transitioning between the two states shown 
in Fig. SJ In actuator mode, the inflating tube could per- 
form work over long axial distances (by pushing a load 
along the axis) or over short transverse distances (lim- 
ited by the diameter of the inflated state) . Axial motions 
could couple to either liquids inside the tube or solids that 
are attached to the tube exterior. Reversing the trans- 
duction, the system could also convert mechanical motion 
into charge, similar to a piezoelectric sensor. Since the 
collapsed state has the smaller capacitance, a charge Aq 
will leave the tube when an external compressive load 
collapses the tube. For example, approximately 200 eV 
of work will collapse a 100 nm length of a (40,40) tube 
held at U = 15 V (for k s = 0.1), producing a charge 
signal of Aq ~ 30e, (Fig. ^ well within the range of 
modern measurement techniques 21-23 . Unlike piezoelec- 
tric crystals whose fractional capacity to elongate Al/l 
is at most on the order of 10~ 3 , bistable nanotube de- 
vices - for axial transport - can have a range of motion 
comparable to the length of the device itself. The effi- 
ciency of this transducer varies considerably depending 
on the operating voltage, but is comparable to that of es- 



FIG. 4: Schematic depicting the evolution of a double-pinned 
tube from a mostly collapsed state at U < (Tdcgcn to a mostly 
inflated state at U > C/dogcn- 



tablished piezoelectric actuators 2 ^ across a broad range 
of operation. 

The time response of these nonlinear nano-mechanical 
systems is governed by the axial speed of the transition 
front between collapsed and inflated states. This front 
moves at v cx V AE®; in our system AE is proportional 
to the voltage deviation from t/degen- For example, a 
ten-volt swing away from £/dcgcn for a (40,40) tube yields 
AE a ~ 2 meV per carbon atom; taking mc as a char- 
acteristic mass, one obtains v ~ \J A~E a /mc ~ 200 m/s. 
A typical device dimension of 0.1 micron then implies 
roughly GHz operating frequencies. 

Devices that operate at the degeneracy point: Of 

particular interest are devices designed to operate as close 
to degeneracy as possible. Such a system is exception- 
ally sensitive to external perturbations that upset the 
balance between inflated and collapsed states. For ex- 
ample, a gas trapped in the interior adds a PV term to 
the free energy, so that equilibrium is obtained when the 
gas pressure balances the energy density associated with 
the volume change between collapsed and inflated states: 
P = AE/AA where A A is a change in the cross sec- 
tion area and AE is measured per unit axial length; (this 
simple analysis neglects fluctuations, which are discussed 
later). For an ideal gas, one obtains the differential sen- 
sitivity of the length x of the inflated portion to changes 
in either atom number n or temperature T: 
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Tuning the voltage to U<i egen , we obtain Ax — > oo. Even a 
single interior gas atom could shift the transition zone by 
a substantial amount, limited by the precision with which 
the system energetics can be tuned. Expanding AE{U) 
around U = f/degen for a (40,40) tube with k s = 0.1, 
we can obtain the response functions to temperature 
{Ax/ AT) or atom number (Ax/ An) as a universal family 
of curves plotted in Fig. [SJ These curves are interpretable 
as either the sensitivity to temperature at fixed number 
of interior gas atoms or the sensitivity to atom number at 
fixed temperature. Other external perturbations can also 
elicit strong responses: for example, dramatic responses 
can be anticipated to surface acoustic wave excitation, a 
highly nonlinear version of the nanotube charge pump- 
ing that was previously examined 2 ^, theoretically and ex- 
perimentally for regular inflated tubes. Adsorbates and 
variations in temperature 2 ^ or mechanical boundary con- 
ditions provide additional experimental handles to tilt the 
delicate balance between configurational states. 

Collapsed-based nano-electromechanical devices could 
also be operated in constant-charge mode rather than 
constant-voltage mode, a distinction similar to that be- 
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FIG. 5: Universal response function to temperature Ax/ AT 
(or atom number Ax /An) for different number of encapsu- 
lated atoms n (or temperature T). The red dashed line marks 

t^degcn • 



tween constant-pressure and constant- volume. At con- 
stant charge the qU term is absent but the inflated state 
continues to be favored at higher charge per unit length. 
Around the degeneracy point the system self-organizes 
into a mixed collapsed/inflated state, with higher charge 
density in the inflated section than in the collapsed one. 
This charge bubble - a giant polaron of sorts - could 
be subject to further experimental manipulation through 
scanned probes or split gates. 



Although huge, tunable responses to perturbations are 
possible in this system, non-uniformities along the tube 
will prevent perfect tuning and hence cut off the divergent 
susceptibility. For example, voltage fluctuations of SU = 
10 mV are characteristic of nanotubes or graphene on 
silica substrates^ (although more uniform environments 
are possible on alternative substrates^ 8 ). An inhomo- 
geneity of SU — 10 mV cuts off dX/dn at ~ 18 nm/atom 
at T = 400 K. Thermal fluctuations should be particu- 
lar important in this nanoscale, one-dimensional system. 
A simple one-degree-of-freedom model (which assumes a 
constant shape to the transition region) writes the en- 
ergy of the system as AEx where x is the location of 
the collapse-inflation transition. Simple thermodynamics 
then implies that (x 2 ) — (x) 2 = ^ . Since the system 
has only one characteristic length scale, the same univer- 
sal divergence results, with fluctuations in the location of 
the collapse/inflation transition of the same size as the 
increment in location induced by a single additional in- 
terior gas atom (Eqn. [3]o) . This distance can be many 
nanometers and could be exploited in fluctuation-based 
device modalities similar to those of biological molecular 
motors. 
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